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I. INTRODUCTION 



Loop quantum cosmology (LQC) is a nonperturbative theory of quantum cosmology 
motivated by loop quantum gravity. Homogeneous and isotropic space-times have been 
studied in considerable detail in LQC now and there is a growing interest in the dynamics 
of linear perturbations as some effects specific to LQC could potentially be detected in 
the cosmic microwave background or in primordial gravitational waves. This possibility 
is especially important as cosmological observations are likely to provide the best testing 
ground for theories of quantum gravity. 

In order to study perturbations around a homogeneous background, it is of course neces- 
sary to understand the background first. This is under control as the LQC of homogeneous 
and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) space-times with a massless 
scalar field have been examined in the flat closed [El |6| and open cases 0, with a 

nonzero cosmological constant 10] and also with a Maxwell field rather than the massless 
scalar field [TTj] , in each case using the so-called "improved dynamics" quantization. Numer- 
ical and analytical studies have shown that the trajectory of sharply peaked states closely 
follows the classical trajectory predicted by general relativity so long as the curvature is far 
from the Planck scale. However, deviations from general relativity appear in a contracting 
universe as the curvature nears the Planck scale, at which point gravity becomes repulsive 
and causes a bounce. In this way, LQC provides a quantum bridge between a classical 
contracting universe at early times and a classical expanding universe at late times. 

In addition, the dynamics of sharply peaked states are extremely well approximated by 
what are called the effective dynamics. The effective equations provide a quantum-corrected 
Hamiltonian constraint that can be treated classically and the equations of motion include 
the modifications due to quantum gravity effects. For example, for the flat FLRW space-time 
the Friedmann equation is modified to 
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where the critical energy density is p c = 0.41ppi [12j. The effective equations have been stud- 
ied for all of the isotropic space-times considered so far in LQC, and in each case they have 
been shown to provide an excellent approximation to the full quantum dynamics of sharply 
peaked states, even in the presence of spatial curvature or with a nonzero cosmological 
constant in which case the model is not exactly soluble. 

Going beyond homogeneity and isotropy, the anisotropic Bianchi models [l~3Tl5| and the 
inhomogeneous Gowdy model 3-18] have also been studied and, although they have not 
been studied in as much detail as the FLRW models, it has been shown that the big-bang 
and big-crunch singularities are resolved in these models as well. The loop quantizations of 
the FLRW, Bianchi and Gowdy space-times have recently been reviewed in some detail in 
Refs. (lJij. 

As it is relatively simple to study quantum gravity effects in the effective theory, it was 
in this setting that linear perturbations were first examined in LQC. The two main types 
of quantum geometry effects that appear in LQC (and loop quantum gravity), inverse triad 
and holonomy effects, have been studied separately in the effective treatments of linear 
perturbations. Inverse volume effects have been studied for tensor 22| |23| , vector 24J and 
scalar perturbations 25|, |26|; and holonomy corrections have been determined for tensor 



22l . 23], vector [27| and scalar perturbations as well. 
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These two different types of corrections are important in different settings: holonomy 
corrections are relevant when the space-time curvature nears the Planck scale, while inverse 
triad corrections become important when a physical length scale becomes comparable to 
the Planck length £p\. Some recent results suggest that the relevant physical length scale 
is the wavelength of the inhomogeneous modes being studied, in which case inverse triad 
effects vary from mode to mode, depending on the wavelength of the mode j28|. (For a 
different point of view, see [31] where it is instead suggested that it is the distance between 
neighbouring nodes in the fundamental spin- network that is the relevant length scale.) 

Although work studying linear perturbations around FLRW space-times in the effective 
theory has been very useful and has led to several interesting phenomenological investigations 
(see e.g. [31-34||), it is important to verify that there does indeed exist a quantum theory 
whose dynamics, for sharply-peaked states, is well approximated by these effective dynamics. 
After all, the quantum theory should be the fundamental theory, while the effective dynamics 
are merely a (very) useful approximation. In this paper, we will address this issue by 
constructing a quantum theory whose effective dynamics are precisely those found in previous 
studies of holonomy-corrected equations for scalar perturbations [28-3(j. The quantum 



theory will also allow for investigations into the dynamics of states that are not sharply 
peaked, something that is impossible in the effective theory. 

Three methods have been proposed so far in order to model small perturbations around an 
FLRW background. The first is to use a lattice where each cell is approximated by an FLRW 
space-time. As the parameters of the FLRW model vary from one cell to another, there are 
inhomogeneities at scales larger than the volume of each cell. If the parameters only vary 
slightly from one cell to the next, then the perturbations are small and therefore this setup 
describes small perturbations on an FLRW background. This is similar to the separate 
universes approach [35|, 361 and was first suggested in the context of LQC in [37| where 
the kinematics were studied, but neither the Hamiltonian nor the scalar or diffeomorphism 
constraints were defined in this work. The second method is to treat the small perturbations 
as a test quantum field on the LQC FLRW background, much as was done for a scalar field 
in 38]; a brief presentation of this programme can be found in 39]. Finally, the third 
method, similar to the second, is to perform a hybrid quantization where the perturbations 
are quantized a la Fock on an LQC background. This technique was pioneered for the 
study of the the inhomogeneous Gowdy models [l6-18| and has recently been used to model 
perturbations in inflation 40 



In this paper, we will follow the first method. We will build a cubical lattice where each 
cell is taken to be homogeneous and isotropic, and there will be some appropriate interactions 
between neighbouring cells in order to obtain the correct physics. The parameters of the 
homogeneous space-times may only vary slightly from cell to cell, in order to ensure that 
one can meaningfully speak of an average background, and small perturbations around it. 
This will require some restrictions on the initial conditions in order to be sure that we are 
in the regime of small perturbations. The kinematics will be quite similar those defined in 



37], although they will be presented in a slightly different way due to some recent progress 
in homogeneous LQC [H 0] . In this paper, we will only consider scalar perturbations as this 
will significantly simplify the model. In Sec. |Vj we will briefly comment upon a possible 
generalization which would allow vector and tensor modes to be treated in a similar lattice 
LQC treatment. 

The main limitation of the lattice LQC model that we shall introduce in this paper is that 
it cannot be used to study perturbations with a wavelength A shorter than £p\. The reason 



4 



for this is the following: in order to have a lattice that describes perturbations up to a given 
wavelength, an appropriate number of cells is necessary. However, for the homogeneous 
LQC model in each cell to be trusted, the physical volume of each cell should remain much 
greater than the Planck volume. This means that it is only possible to describe perturbations 
whose wavelengths remain much greater than £p\ at all times and therefore lattice LQC 
cannot be used to describe cosmological models such as inflation where the wavelength of 
the inhomogeneous modes seen today becomes smaller than £p\ at early times. Nonetheless, 
lattice LQC can be used in many other cosmological scenarios, including the ekpyrotic and 
matter bounce models, where the wavelengths of the relevant modes remain larger than the 
Planck length at all times. Note that an important consequence of this is that inverse triad 
corrections shall always be negligible in lattice LQC as they only become important when 
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The outline of the paper is as follows. In Sec. [TTJ the classical perturbation theory for 
scalar modes (in the appropriate variables for LQC) is recalled from j28|, and we present 
a discretization of the theory. In Sec. IHH the restrictions necessary for wave functions 
to represent small perturbation states are described, and then it is shown that for these 
states the scalar, diffeomorphism and master constraint operators weakly commute with the 
Hamiltonian. The effective theory is obtained in Sec. IIVI and we conclude with a discussion 
in Sec. [V] 

In this paper, we set c = 1 but leave G and h explicit so that it will be easy to see 
whether contributions come from gravitational effects, quantum effects or a combination of 
the two. The Planck length is defined as £p\ = \J Gh. Finally, as the perturbations are 
linear, all terms that are second order (or higher) in the perturbations can be dropped and 
therefore all equations are understood to hold up to first order in perturbations (with the 
sole exception of the Hamiltonian where second order terms are relevant as they generate 
the dynamics of the linear perturbations). 



II. COSMOLOGICAL PERTURBATIONS ON A LATTICE 

In order to study small fluctuations around the flat FLRW cosmological background (with 
topology T 3 ), one allows small departures from homogeneity. A particularly nice coordinate 
choice for scalar perturbations is the longitudinal gauge in which case the metric is given by 
(assuming zero anisotropic stress) 

ds 2 = -(1 + 2^)dt 2 + a 2 (l - 2^)df 2 , (2.1) 

where a(t) is the scale factor and depends only on time while ip(x, t) encodes the fluctuations 
away from the mean scale factor and varies both with time and position. We have chosen 
the line element so that the volume of the 3-torus, with respect to the background metric 
given by ds 2 = dx 2 , is 1. 

The longitudinal gauge is a very useful choice as it can be viewed as the standard FLRW 
metric, though with a scale factor that varies slightly with position. This property leads to 
drastic simplifications on the lattice where each cell is approximated as being homogeneous. 
In the longitudinal gauge, the homogeneous metric in each cell is precisely the FLRW metric 
and therefore, the loop quantization of each cell will be relatively straightforward as the LQC 
of flat FLRW space-times is well known [3, 0]. It is this major simplification that motivates 
the choice of the longitudinal gauge here. In other gauges for scalar perturbations (and for 
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vector and tensor modes no matter the gauge), the metric is necessarily nondiagonal and 
then the well understood results of isotropic LQC cannot be applied in a lattice LQC setting. 

In the first part of this section we will briefly present a canonical treatment of scalar cos- 
mological perturbations in the longitudinal gauge, closely following the results of [28j] . As we 
are only considering linear perturbations, the scalar, diffeomorphism and Gauss constraints 
can be truncated to first order. (In fact, the Gauss constraint vanishes to first order due 
to the choice of the longitudinal gauge and it can therefore be ignored.) However, in the 
Hamiltonian, terms that are second order in the perturbations generate the dynamics of first 
order perturbations and therefore second order terms are relevant for the Hamiltonian. Note 
that the Hamiltonian vanishes to first order, but not to second order. See (41, 42] for further 
discussions regarding these issues. The main ingredients introduced here will therefore be 
the Hamiltonian, to second order in perturbation theory, and the scalar and diffeomorphism 
constraints, to first order in perturbation theory. This will be done using variables that will 
be convenient later for the nonperturbative loop quantization of the FLRW space-times in 
each cell of the lattice. This review will be brief, for further details concerning the canonical 



treatment of cosmological perturbations see [28|, |4l|, |42 



In the second part of this section, we will discretize the Hamiltonian as well as the scalar 
and diffeomorphism constraints on a cubical lattice. We will see that the Hamiltonian 
naturally breaks into a homogeneous term and an interaction term; as we shall see, this split 
will allow us to perform a loop quantization on a lattice in Sec. IIHI 



A. Perturbations in Tetrad Variables 

In order to study the perturbations in the loop gravity framework, it is necessary to use 
the Ashtekar connection and densitized triads as our basic variables. Since the metric is 
diagonal, we can parametrize the densitized triads E% = \Z\q\ef by 

E? = p^e*}, where p = a 2 (l - 2ip), (2.2) 

where p is a function of position and time. We will drop the arguments for the remainder of 
Sec. Ill Al in order to simplify the notation, except in the definition of the Poisson brackets 
where they are essential. We assume that J M ip = in order to ensure that a is the true 
background. Therefore, f M p = a 2 - 

The Ashtekar connection is given by A l a = T l a + jK*, where 7 is the Barbero-Immirzi 
parameter, T l a is the spin-connection and K % a = K a be bl is related to the extrinsic curvature. 
A very useful propery of the Ashtekar connection in this context is that its diagonal terms 
solely come from K l a , while its off-diagonal terms solely come from r* . As the densitized 
triads are diagonal, only the diagonal part of the Ashtekar connection will appear in the 
induced symplectic structure and therefore it is convenient to parametrize the Ashtekar 
connection by 

A x = ct 1 - (dgi/j) t 2 + {dyijj) r 3 , 

A y = {d z ijj) n + cT 2 - (d^) r 3 , (2.3) 
A z = -(dytfj) n + (d x ip) t 2 + ct 3 , 

where A a = A l a Ti and the r, are a basis of the Lie algebra of SU(2) such that TjTj = 
\(-ij k Tk — and we have solved for the spin-connection in terms of the perturbation ip. 
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Following this definition, we find that the induced symplectic structure on our phase space 
gives the following nonzero Poisson bracket: 



{c{x),p(y}} 
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(2.4) 



We will consider the case where the matter field is given by a massless scalar field in which 
case the fundamental Poisson bracket for the matter sector is given by 



{<p{x),n v (y)} = S 3 (x-y). 



(2.5) 



In order to obtain the Hamiltonian constraint, one must first derive the scalar, diffeomor- 
phism and Gauss constraints. For linear perturbations around the flat FLRW space-time in 
the longitudinal gauge, we find that the Gauss constraint is automatically satisfied, while 
the scalar constraint to second order in perturbation theory is 
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2p 3 / 2 g 2 V V 



0. 



(2.6) 



where V is the gradient and A denotes the Laplacian. The diffeomorphism constraint to 
first order in perturbation theory^ is given by 



d a c + c d a ijj 



+ TT^daif RJ 0. 



(2.7) 



Finally, since N = 1 + ip and N a = in the longitudinal gauge, the Hamiltonian constraint 
is given by 



C 



M 



NU + N a U a + A% 
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(2.8) 



See [28j for a derivation of the results presented above. 

We will now introduce a new pair of variables which are convenient for the loop quanti 
zation, 

sgn(p)|p| 3 / 2 a 3 (l-3ip) 



27T74, ' 



and 



which are canonically conjugate, 



p (27T 7 4l) 1/3 |H 1/3: 



{b(x),u(y}} = -5 3 (x-y). 
In these variables, the Hamiltonian constraint is given by 



(2.9) 
(2.10) 

(2.11) 
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(2.12) 



It is only necessary to calculate the diffeomorphism constraint to first order as the shift vector N a is 
necessarily small and as the diffeomorphism constraint appears in the Hamiltonian as N a/ H a , terms that 
are first order in the diffeomorphism constraint will provide second order contributions to the Hamiltonian. 
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where C]^ om is the ultralocal (or "homogeneous" ) portion of the scalar constraint which does 
not include any derivatives and is given by 
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and the interaction terms are 
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where on the first line we introduced 



qv, 



(2.14) 



(2.15) 



M 



and used the relation (1 + ■p)is l /' i = (1 + '0)P 1 / 3 (1 — ip) = is 1 / 3 which holds to first order in 
perturbation theory. Also, the first term on the first line vanishes as j M V 2 ip = (note that 
v can be pulled out of the integral). In the Hamiltonian constraint, terms that are first and 
second order in the perturbations are relevant for the dynamics, while higher order terms 
can be ignored. 

Note that the homogeneous term in the scalar constraint here has the exact form of the 
scalar constraint of the flat homogeneous and isotropic FLRW space-time obtained in (2| 
(with the difference that now b and v are position-dependent.) We will take advantage of 
this in the quantization procedure in the following section. 

This Hamiltonian, expressed in terms of the basic variables b, is, ip, ir^, provides the dy- 
namics. It follows that all of the occurences of ip in Ch must be replaced by 



iP 
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Note that this means that the nonlocal quantity u will appear in the Hamiltonian constraint 
and therefore, in principle, the equations of motion could be nonlocal. However, one can 
check that all of the nonlocal terms vanish and that the resulting equations of motion are 
local, as expected. 

The equations of motion are obtained from the Poisson bracket via 



6{x) = {0(x),C H }, 



(2-17) 



and the equations of motion for the basic variables, truncated to first order in perturbation 
theory and only considering positive values of v for simplicity, are given by 



1 + Ip 7T„ 
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(2.18) 

(2.19) 
(2.20) 

(2.21) 
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In order to obtain the last equation of motion, the relations H — and f M V 2 vp = are 
needed. 

Finally, the initial data is constrained by H and 7i a , truncated to first order in the 
perturbations: 

«=«W + < -P^»A^0, (2,2) 
47 4^7^11/ 1 47rG 



fta = y <9 a 6 + n ^ ~ - ( 2 - 23 ) 



A nice feature of these variables is the particularly simple form of the diffeomorphism con- 
straint. One can check that the scalar and diffeomorphism constraints weakly commute with 
the Hamiltonian, showing that the dynamics preserves the contraint surface. 

This concludes the presentation of the ingredients necessary for a canonical treatment of 
cosmological perturbations in loop variables. 



B. Discretization 

To date, it is only known how to quantize homogeneous models in LQC (aside from 
hybrid quantization schemes). Therefore, it seems that the simplest way to incorporate 
inhomogeneities in LQC is to introduce a lattice where each cell in the lattice is taken to 
be homogeneous but the gravitational field and the scalar field are allowed to vary from cell 
to cell, thus allowing inhomogeneous degrees of freedom at scales larger than the size of the 
cells. By following this procedure, we can use the techniques developed for homogeneous 
models for LQC in order to quantize each cell and thus obtain a quantum theory which 
allows large-scale inhomogeneities. 

For simplicity, we decompose the space into a cubical lattice with Z 3 cells and we label 
each cell by a vector z = (zi,z 2 ,z 3 ) where the z a are integers running from 1 to Z. We 
also introduce the unit cell displacement vectors z a where, e.g., z,\ = (1, 0, 0). Note that the 
distance between the centres of two neighbouring cells with respect to the fiducial metric is 
1/Z and therefore the norm of z a with respect to the fiducial metric is \ jZ. 

The Hamiltonian formalism introduced in the previous section can be defined on the 
lattice in a relatively straightforward procedure. First, the induced symplectic structure has 
changed to 

2Z 3 

{bfr),v&)} = -fr-h lA , (2.24) 

{ V (z l ),n v (z 2 )} = Z 3 5^ 2 , (2.25) 

where 8s u z 2 is the three-dimensional Kronecker delta. It would be possible to rescale v — > Z 3 v 
and ix v Z 3 iT !p in order to remove the presence of the Z 3 in the Poisson brackets, but we 
will not do this in order to explicitly show the presence of the lattice. 

Now, in order to implement the scalar and diffeomorphism constraints on the lattice, the 
derivatives must be approximated by differences. There is a choice to be made concerning 
the exact form of the differences used and we will choose 

9o/(f) L-f (f(z + 2z a )-f(z-z a )) = (V d )«/(z), (2.26) 
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(/(* + 2z a ) - f(z- zSj (/(£+ z a ) - f(z- 2z a ] 



a=l 



(2.27) 



a=l 

- 2f(z- z a ) + f{z- 2z a ) + f(z-4z a j) = A d f(z). 



(2.28) 



This choice is a good one because the discretization of (V/) 2 is such that there are no 
f(z) 2 terms which can cause anomalies in the quantum theory (instead there are terms 
like f(z + z a )f(z — z a ) which do not cause anomalies). Note that simpler discretizations 
are possible classically, but become problematic in the quantum theory. The discretization 
Aif(z) has been chosen as it is the one which is compatible with (V/) d , 



^(v/($) 2 d = -2A d /(^ 



(2.29) 



and similarly (Vd) a / is compatible with A d / as 
Z 3 / 

Ad ^ f ) = e S (( V dU0?+ 24) + (V d ) a /(f) - (V d )J(z- z„) - (V d )a/(^ - 34) 

a=l 
1 3 

= ^((Vd) fl [(V d )J](il + (Vd) a [(V d ) a /](£-2l a )). (2.30) 



a=l 



Note also that (V/) d can be written in terms of (Vd) /- 

As the discrete difference operators approximate derivatives, the product and polynomial 
rules should hold for the difference operators, 



(Vd)a[/^]=/(V d ) a [^]+p(V d ) [/], 

(Vd).[T]^nr l (Vd).[/]. 



(2.31) 
(2.32) 



Of course, these relations are not exact as we are working on a lattice, but they are valid 
within the approximation scheme we are using; they can be derived if the discrete derivative 
(I2.26P is assumed to be small, as it should be in the case of linear perturbations on a 
homogeneous background. This derivation is done explicitly in the quantum setting in Sec. 

mm 

Finally, the discrete versions of the relations (I2.15P and (I2.16p . 



-T 

Z 3 ^ 



v[z\ 



and ip(z) 



3 3u ' 



(2.33) 



will be needed, as well as the condition that ^gipiz) = which ensures that v captures the 
homogeneous mode of v. 
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It is now possible to write the Hamiltonian constraint for the lattice. The homogeneous 
part of the Hamiltonian becomes 



n horn 



-Y 

73 



4 v(z) 



.^.\u(z)\b(z) 2 + 
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and the interaction terms are 

(2 7 r 7 £| 1 |K^)l) 1/3 
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Similarly, the scalar constraint on the lattice is given by 



U{z) 



3h 

47 



+ 



(2 ^' )V3 AMD 



and the diffeomorphism constraint is 

tiv(z) 



n a (z) 



47r 7 4 1 |z/(z)| 12irG\v(z)\ 2 / 3 



(V d )«&(*)+7T v (;s)(Vd)a¥>(z)- 



(2.36) 



(2.37) 



The equations of motion can be obtained from Ch = C H om + Off and it is easy to check 
that they are the discretized equivalent of Eqs. (12. 181) — (I2.21j) . 

It is also possible to determine how the scalar and diffeomorphism constraints evolve with 
time on the lattice: 



U(z) 



1 



h 



(2^74^(^)2/3 

1 

2(2vt74iK^)) 2/3 ^ 



u(z)A A b(z) + TT^zjA^z) 



[(Vd) a n a (z) + (V d ) a H a (z - 2z a ) 



(2.38) 



U a (z) = --±—>H{z){V d ) a v{z) w 0, (2.39) 
3z/(z) 

where we have used the relation (Vd)aV' — — [(Vd) a ^]/3z/ (note that analogous relations exist 
for each of the discrete derivative operators), and the product and polynomial rules ( I2.3ip 
and f l2T32|) . 

This shows that the scalar and diffeomorphism constraints are preserved by Ch, up to 
some small errors introduced by terms that are second order or higher in the perturbations 
and can therefore be ignored. 

The most important result of this section is that the Hamiltonian can be seen as a flat 
FLRW space-time in each cell (although the lapse has an unusual form), with b,i/,(p,7T v 
varying from cell to cell, plus interactions between cells. This will allow us to quantize this 
model in the next section by performing the standard LQC quantization for FLRW space- 
times in each cell, and then work with some interaction terms in the Hamiltonian that are 
relatively easy to handle. 
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III. LATTICE LOOP QUANTUM COSMOLOGY 

In this section, we will quantize the discretized model of linear perturbations on a flat 
FLRW background which was presented in the previous section. Since in the discrete version 
of the theory the metric in each cell is FLRW, it will be possible to quantize the discretized 
theory cell by cell using methods that were first developed for isotropic spaces in LQC 

However, there will be one important difference: in standard LQC, the scalar constraint 
vanishes exactly and the Dirac quantization procedure for constrained systems is imple- 
mented exactly. This is not possible in lattice LQC as even in the classical theory the 
scalar and diffeomorphism constraints only vanish to first order, as does the Hamiltonian. 
Because of this, in the quantum theory it is only possible to demand that the constraints 
approximately annihilate physical states. To be precise, the amplitude A 2 of the linear per- 
turbations of a wave function \l> is determined by demanding that the norm of any relevant 
linear perturbation operator acting on ^ is (at most) of the order of A 2 <C 1. Then, for the 
constraint C to hold to first order of perturbation theory, the norm of C^/ must be of the 
order of A A . 

Also, in lattice LQC there is more than one constraint and therefore the constraint algebra 
is not trivial: it is necessary to check that the scalar and diffeomorphism constraints commute 
weakly with the Hamiltonian, to linear order (i.e., the norm of the terms that are not weakly 
zero acting on the wave function must be of the order of A 4 or smaller). Although there 
are some nontrivial checks to be performed in order to ensure that the quantum dynamics 
preserve the constraints, we do not have access to a full constraint algebra here. This is 
because we are working in the longitudinal gauge and the lapse and shift have been fixed in 
order to preserve the gauge. Therefore, as one cannot freely choose any lapse or shift, the 
full constraint algebra is not available. Nonetheless, the commutators between the scalar 
and diffeomorphism constraints with the Hamiltonian are an important check and show that 
lattice LQC is anomaly free. 

In this section, we will begin by defining the kinematical Hilbert space for the lattice, 
first cell by cell and then for the entire lattice. We will next discuss the conditions states 
must satisfy in order to represent small perturbations, before constructing the Hamilto- 
nian constraint operator which will provide the dynamics and define the physical Hilbert 
space. Finally, we will show that the scalar and diffeomorphism constraint operators weakly 
commute with the Hamiltonian. 

A. The Kinematical Hilbert Space 

The kinematical Hilbert space of lattice LQC is the tensor product of the kinematical 
Hilbert spaces for each cell in the lattice. In turn, the kinematical Hilbert space of each cell 
is the tensor product of the gravitational and matter kinematical Hilbert spaces. 

We shall briefly recall the kinematical Hilbert space for the gravitational sector in ho- 
mogeneous and isotropic LQC for the sake of completeness. The kinematical Hilbert space 
is usually presented in the p representation where the operators corresponding to fluxes of 
the densitized triads and holonomies of the connection are defined. However, in the loop 
quantization, it turns out that the holonomies appearing in the Hamiltonian are those with 
a length of ft = X/\/]p\ which depends on the value of p. (Here A is the square root of 
the minimal nonzero eigenvalue of the area operator in loop quantum gravity and has di- 
mensions of length.) Because of this nontrivial dependence on p, it is useful to work in 
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a different representation where the fundamental operators correspond to v and complex 
exponentials of b. For more details of this, see e.g. [2|. It is the kinematical Hilbert space in 
this representation that will be presented here. 

In the \u) representation, the gravitational kinematical Hilbert space (for one of the cells 
in the lattice) is spanned by wave functions ^(v) that are a countable linear combination of 
the basis states, 

\V) = ^W(u)\u), where ^ |^(z/)| 2 < oo. (3.1) 

V V 

Note that this is a discrete sum, not an integral. The inner product is given by 

(p'\u) = 6^, (3.2) 

where 8 v > tU is a Kronecker delta rather than a Dirac delta. 

As the orientation of the triads (right- or left-handed) is encoded in the sign of u, the 
parity operator II which flips the orientation of the physical triads acts as 

f[#(i/) = #(-*/). (3.3) 



Following loop quantum gravity, the operators which are well-defined in LQC are 
holonomies and areas. Due to the homogeneous and isotropic nature of each cell, it suffices 
to consider holonomies which are parallel to the sides of the cell, and the physical areas of 
the sides of the cell. However, as mentioned above, since the holonomies appearing in the 
LQC Hamiltonian are of the form of exp[±z(A/ \Zjp[)c] [3, 0], it is useful to change variables 
to b and v. Then the first fundamental operator that is to be defined is v which acts by 
multiplication, 

v\v) = v\v), (3.4) 

and is related to the volume (of a cell), 

V\v) = 27r 7 4 1 z>/Z 3 |z/). (3.5) 
The other fundamental operator is the holonomy operator that acts by shifting v, 

exp(ii&)|z/) = \v- 2Z 3 £). (3.6) 

The kinematical Hilbert space for the wave functions *&(<p) describing the matter sector 
is the space of square-integrable functions L 2 (R, dip). The basic operators are defined as 
usual: (p acts by multiplication while fr^ = —ihd/d<p acts by differentiation. 

The complete kinematical Hilbert space of one cell is given by the tensor product of the 
gravitational and matter Hilbert spaces, while the kinematical Hilbert space of the entire 
space is the tensor product of the Hilbert spaces of each cell. 

As the physics is invariant under a change of the overall orientation of the triads, the 
complete state should be invariant under the parity transformation 

n t ot*(^l, V2, ■ ■ ■ , = *(-fl, -U2, • • • , -Vza). (3.7) 

Note however that a typical state will not be invariant under the action of a single Tl(z) as 
this only reverses the orientation in one cell. 
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Finally, it is useful to define the operator corresponding to the mean value of v(zi). The 
operator v is simply given by 

1 ZZ 

r}^(u l ,u 2 ,...,u z3 ) = — ^v i ^{v l ,v 2 ,...,v z s). (3.8) 

8=1 

Similar "barred" operators can be defined for other operators by averaging over all of the 
cells in the lattice, but the v operator is the only "barred" operator that appears in the 
Hamiltonian. 

B. Restrictions on the States 

Since we are interested in studying small perturbations, it is necessary to restrict our 
attention to states where the wave function in each cell only varies slightly from one cell to 
another. We will impose two conditions in order to ensure that the wave functions describe 
small perturbations away from homogeneity. 

The first condition is that the parity cannot vary from one cell to another. There- 
fore, the wave functions may only have support on the two configurations where the 
sign of the arguments of the wave function are all equal, i.e., ^(hl, h|, • • • , \^z 3 \) an d 
-\u 2 \, • • • , -|^s|). It follows that@ 

*(..., N,..., -|z/,-|,...) = 0, ViJ. (3.9) 

Since the wave functions only have support on the sections where the Vi are either all positive 
or all negative, and these two sections are related by the parity transformation 

*(-N, -hi, • • • , -M) = n^*(M, hi, • • • , hH), (3.io) 

which must be conserved by the dynamics of the system, we can restrict our attention to 
the "positive" portion of the wave function where all of the arguments satisfy the condition 
Vi > 0. This is what we will do for the remainder of the paper. 

The second condition is that the perturbations in both the gravitational and matter 
fields must be small and therefore the difference between an operator acting on different 
cells should be small as this is, by construction, a measure of the perturbations around the 
homogeneous background. To be more precise, we ask that the norm of a wave function 
acted on by some difference operator A(z\) — A(z 2 ) be much smaller than the same wave 
function acted on by the operator A(z 3 ), for all zi,z 2 ,z 3 . We can write this condition on 
the wave function as 

(iO?0 - A(z 2 yj v < A(z 3 ) (3.ii) 

where '-C' indicates that the norm of the resulting wave function on the left hand side is 
much smaller than on right hand side. Similarly, we also have 

(A{z x ) - A{z 2 )) (B{zz) - B{zS) * « (A(z 5 ) - A(,? 6 )) (3.12) 



2 This can only be done if the evolution generated by the Hamiltonian operator does not couple the positive 
and negative arguments of the wave function. As we shall see, the Hamiltonian satisfies this property. 
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This last relation holds regardless of the ordering of the two difference operators on the 
left hand side. An important point is that these relations must hold whether the A and B 
operators act on the gravitational or matter parts of the wave function. 

As we are only interested in linear perturbations, we will drop all terms that are second 
order in terms of difference operators, i.e. terms like the one on the left hand side of the Eq. 
(I3.12p above, except in the Hamiltonian constraint operator where they generate nontrivial 
dynamics for the linear perturbations. 

This last condition implies several other useful (approximate) relations that hold so long 
as we work with linear perturbations. The first one is simply 

A{z x ) (B(z 3 ) - B(z A )) * = (a{z x ) - A(z 2 )^j (b(zs) - B(zS) * 

+ A{z 2 ) (b{z 3 ) - B(z 4 j) tf, 
= A(z 2 ) (B(z 3 ) - B(z A )) (3.13) 

where we see that no matter which cell the operator A acts on, the result will be the same 
up to the first order in perturbations. An important point is that this relation holds whether 
A is on the right or the left side of the difference operator. Note also that the prefactor to 
Eq. (I3.13P could be chosen to be the barred "average" operator A instead of A(z 2 ). 

Another useful relation is the quantum equivalent to the product rule that we introduced 
for the lattice in Eq. f l2UTj) : 

([AB){ Zl ) - {AB){z 2 )) * = (A{ Zl ) - A{z 2 )) B{z 2 )m + A( Zl ) (b( Zi ) - B(z 2 )^ ^ 

(a{ Zi ) - A{z 2 )) B(z 3 )V + A(z 4 ) (b( Zi ) - B(z 2 )) (3.14) 

where on the second line we used the relation (I3.13P twice. Note that this relation follows 
from the condition (I3.12p . it is not a further assumption. This last relation immediately 
implies a relation analogous to the polynomial derivative, 

(A{z 2 ) n - i(^ 3 ) n ) * = nAiz^- 1 [A{z 2 ) - A(zs)) vl/, (3.15) 

which holds for fractional and negative n so long as A is invertible. The existence of these 
relations for difference operators — that are very similar to those for differential operators — 
should not be surprising as we are only interested in linear perturbations and we are dis- 
carding all higher order terms. 

The final relation we will present here is obtained by 

{AB){z x ) (C(z 3 ) - C{zS) * = A( Zl )B(z 2 ) (C(z 3 ) - C(z^ * 

= B{z 2 )A{ Zl ) (C(z 3 ) - C{zS) * 
= {BA){z 2 ) (c{z 3 )-C{z A ))^ 
^{BA){ Zl )(c{z 3 )-C{z A ))^, (3.16) 
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where we have used the relation (13 . 13[) repeatedly and, on the second line, freely commuted 
the A and B operators as they act on different cells (this requires the assumption that both 
of the operators act only on one cell). What this calculation shows is that terms of the form 

[AB}( Zl )(d(z 3 )-C(z,))^, (3.17) 

even though they may not look as though they are second order in perturbations, can in 
fact be neglected. This relation is equivalent to 

[A,B]{z)V < {AB){z)V, (3.18) 

which, although it at first seems to be quite different from ( 13. 12ft . is in fact implied by it. 

A consequence of the relation (13.18p is that by choosing A = e~ l2Xb and B = v — two 
operators that play an important role in the Hamiltonian — we obtain the condition 

4AZV^A6(z)# < e^(z)v{z)m, (3.19) 

which becomes 

4AZ 3 ^ < v{z)¥. (3.20) 

when it is multiplied by the operator e l2Xb (z). This shows that the wave function must have 
support mostly on large v. 

The condition (I3.20p has a clear physical meaning: it means that lattice LQC can only 
be used so long as the physical volume of each cell remains much larger than the Planck 
volume [recall the definition of the volume operator for each cell given in Eq. (13. 51) ]. As 
the wavelength of the shortest inhomogeneous mode that can be studied in a given lattice 
is determined by (twice) the average physical length of each cell, it follows that, in lattice 
LQC, the wavelength of each mode studied in lattice LQC must remain much larger than 
£p\ at all times. This is a strong drawback as it indicates that, for example, lattice LQC 
cannot be used if there is inflation. Nonetheless, lattice LQC can be used in many other 
cosmological settings, including the ekpyrotic scenario or matter bounce models. 

Another important consequence of the wavelength of all modes remaining much larger 
than the Planck scale is that in lattice LQC inverse triad corrections will be completely 



negligible. In [28], it is shown how to incorporate inverse triad corrections in a lattice LQC 
setting and the main result is that all nontrivial corrections depend on the ratio of £p\ to 
the wavelength. Therefore, if the wavelength is always much larger than the Planck length 
(as it must be in lattice LQC), inverse triad effects are completely negligible. 

As an aside, we point out that the condition (13.201) can also be derived in the following 
manneiH. In the construction of the nonlocal curvature operator used to define the Hamil- 
tonian, one takes a holonomy around a square loop of area ~ i\ v In order to ensure the 
straight edges defining this holonomy remains within one cell, the condition that the length 
scale characterizing the cubical cell be larger than tp\ must be satisfied, which is equivalent 



3 There is also another way to understand (|3.20[) . coming from the Hamiltonian (|3.25p defined in Sec. MI Cl 
together with some results regarding the superselection sectors in Sec. MIDI Restricting our attention to 
positive v and taking the superselection sector e = 4 for simplicity (this argument can easily be extended 
to any other superselection sector), the wave function ^(y) can only have support on v — AmXZ 3 , with 
7?i G N. Now, for perturbations to be small, their amplitude must be much smaller than the mean value 
of the background. Due to the superselection sectors, the smallest amplitude of perturbations possible in 
v is 4A.Z 3 . Therefore, for the perturbations to be much smaller than the background, the mean value of 
v must be much larger than 4XZ 3 . 
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to demanding that the volume of the cell be larger than £ Pl , and this is precisely (13.201) . 
Considering the same argument from the point of view of the continuum (rather than the 
discrete theory), condition (I3.20p can be understood as being equivalent to the assumption 
that the connection is varying slowly enough so that gradients of the connection are neg- 
ligible for holonomies of length £p\. This clearly holds for linearized perturbations whose 
wavelengths are larger than £p\. Note however that this may no longer be the case for 
short (< £pi) wavelength modes in which case the gradient of the connection could become 
important and a different quantization procedure may be necessary. 



C. The Hamiltonian and the Constraint Operators 

In the classical theory, the Hamiltonian as well as the scalar and diffeomorphism con- 
straints do not exactly vanish and therefore it is impossible to ask that their operator analogs 
annihilate physical states in the quantum theory. Instead, following the classical theory we 
will impose that the Hamiltonian and the constraints must be zero to leading order in the 
perturbations. This can be done by determining the strength of the perturbations (which 
must be small) by calculating the norm of a wave function \]/ (satisfying the conditions 
in Sec. IIIIBI) acted upon by an operator linear in perturbationsj. Denoting this norm by 
A 2 < 1, the condition that a constraint C annihilates the state \l/ to linear order in the 
perturbations is that 

\\C^\\ 2 = 0(A 4 ). (3.21) 

Now that it is clear what is meant by a constraint annihilating a state to leading order 
in the perturbations, it is possible to define the relevant operators corresponding to the 
Hamiltonian and two constraints, all of which annihilate physical states to first order in 
perturbation theory. 

In order to define the Hamiltonian operator for lattice LQC, we will follow a simple 
prescription: we will quantize the homogeneous terms in the Hamiltonian, given in (I2.34[) 
just as is done in homogeneous LQC [H, y] and the interaction terms (12.351) will be quantized 
in a simple manner determined by their form in the unique consistent formulation of effective 
equations in the longitudinal gauge [28, Hi. We will start with C^ om , which is easy to quantize 
as it is basically the Hamiltonian in flat FLRW models, although the lapse, N = 4/3 — v/?>v, 
has an unusual form. 

In the loop quantization, there are two nontrivial steps that are necessary in order to 
obtain the Hamiltonian operator in a homogeneous and isotropic space (which each cell 
in the lattice is). First, as the basic operators in LQC correspond to holonomies of the 
connection, it is necessary to express the field strength of the Ashtekar connection. This is 
possible as the field strength can be expressed in terms of the holonomy of the connection 
around a small square loop. Due to the discreteness of the area eigenvalue in loop quantum 
gravity, it is posited that it is appropriate to choose the area of the square loop to be the 
smallest area eigenvalue in loop quantum gravity, denoted by A 2 = A\/?>ti^£ vv A careful 



implementation of this procedure 0, 1 |13| shows that we should express the field strength as, 
e.g. 

iV = (2?T7 A t )2/3 (^W Aft), (3.22) 



4 Examples of such an operator could be the third term in the scalar constraint (|3.29j) (the first two terms 
contribute to both the zeroth and first orders) or either of the terms in the diffeomorphism constraint 
(1X30)1 . 
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where we have not yet specified a factor-ordering on the right hand side. 

The second nontrivial step is to define an inverse triad operator. This is necessary as 
the state \u = 0) is normalizable and therefore zero is in the discrete spectrum of the 
v operator. This difficulty can be avoided in loop quantum gravity by using one of the 



Thiemann identities introduced in 43] and a similar approach can be used in LQC Although 
there is an ambiguity in LQC as there are many ways to construct inverse triad operators (a 
natural choice for inverse triad operators in lattice LQC models is proposed in 2fj), all of 
the inverse triad operators considered so far in the LQC literature share two properties: first, 
they annihilate the state \v — 0) and second, their eigenvalues are very close approximations 
to \/u for states \u) where V £p { , or, equivalently in lattice LQC, v 3> £p y Z 3 . Of course, 
there can be important inverse triad effects for small, nonzero v which may potentially have 
observational consequences. However, as in lattice LQC we necessarily have v S> £p\Z 3 [see 
f)3.20p ]. inverse triad effects are negligible and therefore we will define the simplest possible 
inverse triad operator that satisfies the two properties above: 

"5". x f° ifi/=0, . 

T-r k> = S i i x i • 3.23 
\u\ IrI 1 / otherwise. 

All other inverse triad operators of different powers can be obtained by raising the operator 
defined above to the appropriate power. Similarly, it is possible to define an inverse triad 
operator for the mean operator u, 

T fo if|p|* = 0, 

7^* = S 1 t • 3.24 

\u\ otherwise, 

where the first case occurs only if \I/ has no support on ^ / 0. Note that it is important to 
take the absolute value as \& is symmetric in v as can be seen from Eq. f)3.7p . 

This choice for the inverse triad operator is a simple one, and while it amounts to ignoring 
potential effects coming from a "true" inverse triad operator, it turns out that the effect of 
any other inverse triad operator is extremely well approximated by 03. 23!) as lattice LQC 
only provides a valid description of linear cosmological perturbations if the wavelengths of 
these perturbations remains considerably larger than the Planck length in which case the 
observational consequences of inverse volume operators are completely negligible as argued 
above. 

Now that the field strength has been expressed in terms of holonomies of the connection 
and the inverse triad operators defined, it is possible to define the Hamiltonian 

Ch = C^ + (3.25) 

in terms of the homogeneous and interaction terms. Choosing a factor-ordering where the 
lapse is placed on the left and expressing the sin 2 Xb operator in terms of complex exponen- 
tials, 



nhom 




3fr 



(3.26) 



I67A 

h IT d 



47r-fG\u(z)\d^(z) 2 _ ' 
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This is a nice factor-ordering choice as it ensures that, under the action of C]^ om , the zero 
volume states decouple from the nonzero volume states. It also does not couple the z/j > 
and Pi < sectors (for the relevant superselection sector), which is necessary in order to 
preserve the condition (13.91) and also allows us to restrict our attention to one of the sectors, 
say Vi > 0, when we are working with this part of the Hamiltonian operator. Then, the 
sectors where some or all of the v\ < are given by Eqs. ( 13. 71) and ( 13.9H . We will return to 
these points at the beginning of Sec. IIII Dl 

Recent studies of holonomy-corrected effective equations for scalar perturbations in the 
longitudin al g auge indicate that there should only be one modification to the inhomogeneous 



terms [28|, |29|: the second term in Eq. (12 .35j) should be multiplied by cos2A6. This gives 



z 



1 1 

Vviz 



9v(z) 2 V ' v "V d 

u(z)\^, (3.27) 



+ 47rGcos2Xb (V0(z)y 

where we have chosen a symmetric factor ordering in the last term. The difference terms 
are defined in Sec. IIIBj for example 

/ \ 2 Z 2 ^ / \ / \ 

(V0(2)J = -j- J2 Tz -) - *«)) + 2 «) - <P(z- 2««)J. (3.28) 

a=l 

By following the same steps as above, the expression for the scalar constraint operator is 
given by 

H(z) =^^/W)\{e : < Si (z) + e^(z) - 2) VR^| 



" 1 * ^^..T^ l^l, (3,9) 



47r 7 G|z/(i')|^(f) 2 12ttG |z/(z)| 2 / 3 

where the factor-ordering is chosen following C|f m . 

The diffeomorphism constraint operator is slightly more difficult to define. This is because 
an effective analysis of the model shows that the diffeomorphism constraint should not be 
modified, even when the curvature is of the Planck scale j28[[29[. Nonetheless, it is necessary 
to express all occurances of b in ( I2.37P by a complex exponential. Luckily, there is an easy 
way to ensure that the diffeomorphism constraint is not changed in the effective limit which 
is by also multiplying the matter sector of the diffeomorphism constraint by the appropriate 
complex exponential of b, giving 

H a (z) = !^(V d ) a ^(;?) + ^(z)7c v (z)(V d ) a 0(z). (3.30) 

The choice of the factor of 2 in the complex exponential is necessary in order to have an 
anomaly-free constraint algebra in the quantum theory. [See Eq. ( I3.63|) and the calculations 
leading to it for an a posteriori justification of this.] An important property of the diffeo- 
morphism constraint defined above is that, since it annihilates all physical wave functions 
(to leading order in perturbation theory), so does the operator 

W a (z) := e^ b (z)U a (z) = -!^©(V d ) a e"SA6(z) + e^ b (z)ri v (z)(V d ) a 0(z), (3.31) 
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where we have used the relation f !3.14j) twice in order to obtain the right hand side of the 
equation above. Therefore, the diffeomorphism constraint operator can be defined as either 
Eq. (I3.30p or (I3.3ip as they are equivalent. 



Finally, it is necessary to define the master constraint Cm |44| as the scalar and diffeomor- 
phism constraints, evaluated at neighbouring lattice sites, do not commute with themselves 
or each otheiH. This can easily be done with the help of the inverse fiducial metric q ab , 

Cm = Y, (u{z)U{z) + q ab U a (z)U b (z)^ . (3.32) 

z 

Physical states in lattice LQC are those that are annihilated, to linear order in the pertur- 
bations as defined above, by the master constraint. 

The Hamiltonian operator and the scalar, diffeomorphism and master constraint opera- 
tors are enough to define the theory of lattice loop quantum cosmology. In the following 
part we shall study the salient aspects of the dynamics and show that all of the constraints 
weakly commute with the Hamiltonian. In Sec. [IV] we will show that the resulting effective 
equations — for holonomy corrections to scalar perturbations in LQC — are the same as 



those obtained in previous works [28|, [29 



D. The Dynamics of Lattice LQC 

A nice property of the Hamiltonian operator is that it automatically resolves the singu- 
larity. If one starts with a nonsingular state, i.e., a wave function that does not have any 
support on v{z) = for all z, then the wave function will continue to be nonsingular. 

This can easily be seen as it is only the shifted terms that allow the support of the wave 
function to change. However, each of terms where the argument has changed to v{n) ±AXZ 3 
has a prefactor of \v(n) ± 4AZ 3 | (to some positive fractional power) which is clearly equal 
to zero if the new argument of the wave function is zero. Because of this, a wave function 
with no initial support on singular states necessarily continues to be nonsingular. It is in 
this precise sense that the singularity is resolved by the Hamiltonian. 

Also, just as in the flat FLRW case, there are also superselection sectors for each cell: if 
one starts with a wave function that only has support on a particular superselection sector, 
then the wave function will continue to only have support on this superselection sector 
throughout its evolution with respect to Ch- The superselection sectors for u(n) are denoted 
by e and are given by € + 4mXZ 3 and A—e+AmXZ 3 for e G (0, 2XZ 3 ) or e + 4mXN 3 for e equal 
to 2XZ 3 or AXZ 3 . In both cases, m G Z. The shorter stepsize in the superselection sector 
for e G (0, 2XZ 3 ) is due to the reflection conditions at v = due to the parity symmetry 

tf(M) = *(!/). 

The most interesting superselection sector is the e = 4 one as this is the only superselection 
sector that includes the singular state v = (although the singular states do decouple in 
any case due to the form of the Hamiltonian). This superselection sector is also particularly 
interesting as it separates the v > and v < sectors: it is easy to see that there is no 
interaction between these sectors for e = 4XZ 3 . Since this is the only superselection sector 
that contains the singular state, and it has the nice property that it separates the positive 



5 Note that there is an important difference between this master constraint and the one proposed in 44j : 
the lattice LQC dynamics are not generated by the master constraint, but by a separate operator, Ch- 
Thus, the role of the master constraint is more restricted here. 



20 



and negative v sectors, we will work with this superselection sector in each cell for the 
remainder of the paper. 

Since the positive and negative sectors decouple, it is sufficient to only consider the part 
of the wave functions where their arguments are positive. The negative part of the wave 
function can be obtained by the parity relation (13 .TP . In addition, assuming that we start 
with a nonsingular state, the \v{z) = 0) singular states can be removed from the Hilbert 
space as they decouple under the dynamics. An important ramification of this is that, once 
the zero volume states have been removed from the Hilbert space, 



v(z) 



viz. 



(3.33) 



These simplifications will make it much easier to write out the action of the Hamiltonian. 

Let us now define our notation: we shall denote states by \I> := (pi) where the index 
i is a shorthand for the vector z denoting the cells in the lattice. For simplicity of notation, 
the arguments of the wave function shall be dropped except when one of the arguments in 
the wave function has been shifted. In this case, only the shifted argument will be displayed, 
e.g., ^ + 4AZ 3 ). 

Using the notation described above and restricting our attention to the e = 4 superse- 
lection sector and v > 0, the action of the Hamiltonian operator is Ch^ 
where 



^hom^ + c 



int \T/ 
H w > 



73 



Z 3 







(1- 


v J 



3h 



(y/vf + AXZ 3 v l ^(v i + AXZ 3 ) - 2u& 



+ 



AXZ 3 Vi^(vi 



AXZ 3 



h d 2 
AnjGvi dip: 



(3.34) 



and 



1 

z 3 



z 3 



i=i \ 



+ fa 



( (2tt 7 £ 



2 \V3 



5/3 



72nGv; 
- AXZ 3 ) 1/6 y(vi 



d 



fa + AXZ 3 ) 1/6 ^fai + AXZ 3 



AXZ 3 



(3.35) 



where the difference operator is defined in (I2.27P — recall that the index i is simply shorthand 
for z — and it is understood that 

1 Z 3 

-V = ^ — (3.36) 

It is also possible to explicitly write out the actions of the scalar and diffeomorphism con- 
straints, but that will not be necessary for our purposes here. 



The physical inner product can be obtained by the group-averaging procedure |45H47 
(assuming the Hamiltonian is essentially self- adjoint), 



da(#|e iQC "|$), 



(3.37) 



while a family of Dirac observables is given by relational observables with respect to the 
massless scalar field ip(z) (48 Hoi . 
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Alternatively, it is possible to deparametrize the system with respect to the massless scalar 
field following 5l[ by using the diffeomorphism constraint in order to write the Hamiltonian 



solely in terms of geometric operators and 7t v (z). Then, solving for fr v (z) allows us to use 
the scalar field <p(z) as a time variable for each cell in the lattice. This "emergent" clock 
automatically provides a physical inner product and an explicit definition of the relational 
Dirac observables. We leave this deparametrization as an exercise for the reader. 

Now that the action of the Hamiltonian operator has been written out explicitly, it is 
perhaps easier to see some of its properties that were described above. In particular, it is 
now more obvious that the singular states Vi — decouple from the nonsingular ones under 
the action of Ch, and also that the positive and negative sectors do not interact. Finally, 
the presence of the superselection sector v = 4m\Z 3 , with m 6 N (for v > 0) is clear. 

Finally, in order for the evolution generated by the Hamiltonian operator to be compat- 
ible with the scalar, diffeomorphism and master constraints, it is necessary to check that 
the dynamics are anomaly-free, i.e., that the Poisson bracket of the constraints with the 
Hamiltonian are weakly zero to leading order in perturbation theory. 

In order to calculate the relevant commutators, it will be necessary to use several of the 
relations obtained in Sec. IIIIBI in order to show that the Hamiltonian and the constraints 
commute weakly. In particular, the relations f)3. 13j) . (|3.14|) . ( 13. 15j) and (13.201) will be used 
repeatedly. We will also only consider states that have no support on v = as these states 
can be removed from the Hilbert space due to the form of the Hamiltonian. (Recall that 
in this case the inverse triad operator corresponding to 1/u is the inverse of v.) Finally, we 
will drop the hats on the operators in order to simplify the notation. 

Because of its nonlocal nature, commutators where v appears require some care and 
therefore we will derive the following commutator — which as we shall see will be a contribu- 
tion coming from the lapse — in some detail. The relevant nonlocal commutator to calculate 
is 



^— ' 3v 



viz) ± 4AZ 3 viz) 



3(P±4A) 3v 

y{y) _ K£> 
3(^±4A) 3v 



e Ti2Xb (z) 



E 



e Tt2Xb {z); 



and it is important to be careful as there is one term in the sum whose commutator is 
different from the others. Then the top term on the right hand side can be separated into 
a part which combines with the nonlocal term and a remainder, this gives 



3(^±4A) V ; ^ 3v{v±\X) y ' 

= T^fe^\z)l ±^J2v(y)e^(z)l (3.38) 

y 

An important observation here is that there is a nonlocal term on the right hand side due 
to the nontrivial commutator of v~ l and the shift operator which could be problematic. 
However, as we shall see shortly, although such a nonlocal contribution does indeed appear 
in the commutators between the Hamiltonian and the constraints, it will be multiplied by 
the scalar constraint (modulo a term which sums to zero) and therefore does not pose any 
problems. 
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The commutator between the scalar constraint and the Hamiltonian is given by 



Kg) 

3v 



+ Z 3 S (3 



4 z/(£) 



int(l) 



3// 



+ 



U\y)+U 2 



[n 2 {z),H\y)] + [U'{z),U l 



H 2 (z),C^ 



(3.39) 



where the superscripts indicate the first, second and third terms of each operator as defined 
in f)3.26p . ( I3.27P and ( I3.29p . There are several terms that obviously commute (either exactly 
or to first order in the perturbations) and these terms have been dropped. In addition, it is 
easy to see that the first two terms on the second line cancel. We shall now calculate the 
remaining terms one by one, to first order in perturbation theory as all higher order terms 
can be dropped. 

The first term can easily be computed by using the result (I3.38p . and gives 



73 L^i 



Z 3 



U\z), 



Kg) 

3v 



n\y) + n 2 
%h\ viz 



2 7 A 

^z/(y)sin2A6( 



,,[ " ] sin2Xb(z)^ 



v 



U\z)+U 2 {z)) 



'viz 



V 



2~i\Z 3 v 

y 

The second line of this result can be simplified by noting that 



(n\y) + H 2 



v{y) sin2A6(; 



-U 3 



viz) sin2A&(z)^-^ V H 3 iy) = 



(3.40) 



(3.41) 



v y 

The first equality holds to first order in the perturbations as there is a difference operator 
in H 3 , while the second equality holds as the sum over all cells of a difference operator van- 
ishes (this is the discrete equivalent of Stokes' theorem on a manifold without boundaries). 
Therefore, by adding this particular zero to the commutator, we find that 



-y 

z 3 ^ 



n l (z), 



Kg) 

3v 
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U[Z] sin2Xb(z)^l(n\z) + n 2 (z] 



27A 



^-^,(2/) S in2A6(,)- 



-H(y), (3.42) 



where it is now clear that the nonlocal term is weakly zero. 
The second nontrivial commutator, 
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3v 
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(3.43) 
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splits naturally into two separate terms. The first term will pick out 6 contributions from the 
sum due to the structure of the difference operator (12.281) (that can themselves be expressed 
in terms of Ad). One might be tempted to drop the second term as it is of the form (13.161) . 
but as the prefactor is very large [see (I3.20p and remember the form of TL ], this term is 
relevant. A straightforward calculation shows that the first term gives 

^^P m2M ' (3 - 44) 
which has been rewritten on the right hand side by using the relations (13.131) . f !3.14j) and 

dSlSD. 

The second term can also be calculated rather easily. It is easy to see that the contribution 
from the second term is given by 



12nGZ 3 \v(z) 2 /z \u{z) +4AZ 3 ] 2 /3 



PI) L x \ -i2\b 



e- i2Xb (z)A d v(z) 



+ 12vrGZ 3 ' [u(z) - 4AZ 3 ]V3 ) e W±im, (3-45) 

which can be simplified by Taylor-expanding (z/±4AZ 3 )~ 2 / 3 to second order [all of the higher 
order terms can be neglected due to (13.20)) ] . This gives 

z^(2 7 r 7 4 1 ) 1/3 sin2 A6(z)A d z/(z), (3.46) 



127tG^Xp(z) 2 / 3 



which will cancel part of the contribution from the first term. 
Combining the results (13.441) and ()3.46p . we find that 



There are two commutators that remain, both involving the interaction terms in the 
Hamiltonian. These are both relatively easy to calculate. By using the relations (13. 13)) . 
(13. 14)) and (13 .15)) . it is possible to express these commutators as 



U\z)^ l) ] =^{4^sin2Xb(z)^n%z), (3.48) 
J 27A v 



H 2 (z),C^ ti2) 



1,1 cos2\b(z)ir ip (z)A d (p(z). (3.49) 



(27T 7 4 1 ) 2/3 ^(^) 2/3 

The awkward factor-ordering in ( 13.481) is chosen to coincide with that in ( I3.42j) . 
By combining Eqs. 1)522)1 . 1EQ7|1 . ( BHD and (E3HD, we find that 

y 
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thus showing that the scalar constraint weakly commutes with the Hamiltonian to first order 
in perturbation theory. 

The next step is to calculate the commutator between the diffeomorphism constraint H, a 
and the Hamiltonian [it is easy to see by the first equality in (I3.3ip that if H a commutes 
weakly with the Hamiltonian, so will W a ]' 



[H a (z) : C H ]^H^,CH] (V d ) a e 



- 2X \z) + ^u(z) [(V d ) B e- 



l2Xb (z) 



C H ] 



+ ^(*)((V d )«¥>(2)J [e- i2X \z),C H ]+^{z)[{V d )Mz),C H ]e- t2X \z) 

+ K{z),Ch] ((V d ) aV (z))e^ 2Xb (z). (3.51) 



In order to complete this calculation, we shall compute one by one the commutators of 
<p(z),7r v (z),v(z) and exp(— i2Xb)(z) with the Hamiltonian. The second and third of these 
commutators need only be calculated to leading order due to the presence of a difference 
operator multiplying each of the appearances of these commutators, while the other two 
commutators must be calculated to first order in perturbations. 

The two commutators that are only needed to leading order are both easily calculated, 
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3h 



47A 



^MF)(e- l2Xb (zi 



J2Xb/ 



K(z),Ch] = 0. 



(3.52) 
(3.53) 



Of the two remaining commutators, both of which must be calculated to first order in 
perturbations, [</?, Ch] is the easier calculation and gives 
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(3.54) 



The last commutator, 
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(3.55) 



is a little more complicated and shall be calculated term by term. Note that although 
C# does not exactly commute with the shift operator, it does up to first order in the 
perturbations. 

The first term (1/Z 3 ) ^(1 + if>(y)) [e- i2Xb (z),U\y) + U 2 (y)] gives 
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while the contribution from the second commutator is 



-i2Xb(->\ pmti 1 ) 
\ Z )i U H 



4A 
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■H 3 (?)e- t2Xb (z). 



(3.57) 



Note that the factor-ordering can be changed freely due to relation ( 13 . 16[) as there is a 
difference operator in H 3 . 

Finally, the last commutator gives 
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and it is possible to combine all of these terms which together show that 
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(3.59) 



where we have changed the factor-ordering in (I3.57P by using the relation (13. 16[) and we have 
added zero to the nonlocal term in (I3.58P in the form of 
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(3.60) 



in order to show explicitly that the nonlocal term is zero on the constraint surface. This 
expression vanishes as the sum of a difference operator over the entire lattice is identically 
zero. 

It is now possible to complete the calculation of the commutator, started in Eq. (I3.5ip . 
between the diffeomorphism constraint and the Hamiltonian by using the results (I3.52p . 
f )3.53p . (I3.54p and ( I3.59p . We will not go through the calculation step by step as it is not 



particularly difficult (although a little tedious), but rather present two relations that are 
necessary for this derivation that may not be immediately obvious. The first, which can be 
derived from the relation (j3.14j) . is 

yjv{2) + AXZ 3 (V d ) a ^u(z) + AXZ 3 = ^(V d ) a (u{z) + 4AZ : 



^d)aV(z) 
V^)(V d )aVvl/) 



(3.61) 



and therefore the difference between the term on the left and the bottom term on the right 
is zero. Two other terms that can also be seen to be cancel by using relation (I3.14p are 
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After a long but straightforward calculation where the relations (I3.13p . ( I3.14p . ( 13 . 1 5[) and 
(I3.16P are used repeatedly, as well as the two equalities derived above, the result 
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(3.63) 



is obtained, showing that the diffeomorphism constraint weakly commutes with the Hamil- 
tonian to leading order in perturbation theory. 

These calculations show that in lattice LQC the quantum dynamics generated by the 
Hamiltonian preserve the scalar and diffeomorphism constraints, and therefore the master 
constraint also weakly commutes with the Hamiltonian. This shows that lattice LQC is 
anomaly-free. 

There are two main limitations here. The first is that we are working in a perturbative 
setting and therefore we can only show that lattice LQC is anomaly-free to first order in 
perturbations. Nonetheless, this is the best that can be done when working with linear 
perturbations, as we have been from the very beginning. The other limitation is that the 
physical volume of each cell in the lattice must remain much larger than the Planck volume 
at all times [see condition ( I3.20p ]. The main consequence of this is that short- wavelength 
modes (with respect to Ep{) cannot be modelled in lattice LQC, and therefore inverse triad 
corrections cannot be studied here. One would expect contributions from inverse triad 
corrections to play a potentially important role for short wavelengths, but to understand 
these it will be necessary to use a different approach to study linear perturbations in LQC 
than the lattice model presented in this paper. 



IV. EFFECTIVE EQUATIONS 

The effective theory is obtained by studying the dynamics of states that are sharply 
peaked. The resulting effective equations should provide the leading order quantum cor- 
rections to the classical dynamics. In general, we do not expect the effective equations to 
be valid when quantum effects are at their most important, but surprisingly the effective 
equations provide excellent approximations to the full quantum dynamics of sharply peaked 
states — even when quantum effects are strongest — for FLRW models. 

As in lattice LQC each cell is homogeneous and isotropic, and the interactions between 
nearby cells are small, it is reasonable to hope that the effective theory will be just as good 
an approximation to the fully quantum dynamics of sharply peaked wave functions as in the 
exactly homogeneous case. 

We will first obtain the effective equations for lattice LQC and then the continuum 
limit will be very easy to take, giving the holonomy-corrected effective equations for scalar 
perturbations in LQC. (Recall that inverse triad corrections are not included due to the 
condition f)3.20p which restricts the wavelength of the modes considered in the lattice to 
always remain much greater than the Planck length, in which case inverse triad corrections 
are expected to be negligible.) 
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A. The Effective Theory on the Lattice 



The Hamiltonian for the effective equations is given by the sum of the terms corresponding 
to the homogeneous and interaction parts of the Hamiltonian given in ( 13.261) and ( 13.271) . 
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(4.1) 



except that Ch (and its components) are now treated as classical objects. The scalar and 
diffeomorphism constraints are similarly the classical equivalents to the operator equations 
(I3T29|) and (ETSOjl . 
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(4.2) 



(4.3) 



Note that the e _?2Ab terms cancel in % a and therefore the diffeomorphism constraint in the 
effective theory is the same as the classical discretized diffeomorphism constraint (12.371) . On 
the other hand, the scalar constraint is different in the effective theory. 

The effective dynamics are obtained by following the same procedure as usual with the 
Hamiltonian and constraints given above. Recalling the Poisson brackets on the lattice 
(I2.24p and (12 .25 p . we find that the effective equations are given by 
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where it is understood that 
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(4.4) 
(4.5) 
(4.6) 
(4.7) 

(4.8) 



It is possible to show that these equations of motion conserve the constraints, i.e., the 



Hamiltonian weakly commutes with the scalar and diffeomorphism constraints j28l 129 



B. The Continuum Limit 

The lattice is essential in the quantum theory in order to have available all of the LQC 
tools that have been studied in considerable detail for homogeneous space-times. Further- 
more it is not clear how to obtain, as a limit from lattice LQC, a quantum theory defined 
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on a continuum, or even if such a limit exists. However, the situation is very different in the 
effective theory where the continuum limit is trivial. 
The effective Hamiltonian, in the continuum, becomes 
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Similarly, the scalar and diffeomorphism constraints are 
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(4.11) 



in the continuum. The equations of motion can be obtained from the Hamiltonian ( 14. 9p and 
using the Poisson brackets ( 12. 51) and 112. lip , giving 
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(4.12) 
(4.13) 
(4.14) 
(4.15) 



In the (continuous) effective theory, the scalar and diffeomorphism constraints ( I4.10p 
and (14. lip must be satisfied and then the evolution of the basic variables generated by the 
Hamiltonian ( I4.9p (which weakly commutes with the constraints) is given by Eqs. (14. 12p — 
( I4.15p . This picture is nicely consistent with previous studies of the holonomy-corrected 
effective equations for linear cosmological perturbations as these equations are in exact 
agreement with those obtained in 28, 

Finally, the correct classical limit is easily obtained from these equations by sending 
H — > (recall that there is an h in A, and see [28J for the classical equations of motion for 
the variables b and u). 



V. DISCUSSION 

Lattice LQC is an anomaly-free loop quantization of linear cosmological perturbations. It 
includes the same quantum geometry effects that cause the bounce in homogeneous models 
and it is easy to see that the singular states decouple from nonsingular ones under the 
action of the Hamiltonian, thus resolving the classical singularity. In addition, the effective 
dynamics of lattice LQC are the same as those derived for holonomy corrections in Refs. 
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28l . [29J, and the correct classical limit is obtained when we take h — >■ in the effective 
equations. 

We expect the effective equations to provide an excellent approximation to the full dy- 
namics of sharply peaked states, just as in the homogeneous case, because in lattice LQC 
each cell is taken to be a flat FLRW model which interacts only weakly with its neighbours. 
Therefore, it seems reasonable to expect that the effective equations will continue to be very 
good approximations to the quantum dynamics of sharply-peaked states. Nonetheless, this 
hypothesis should be tested; this would probably be easiest to do numerically. Numerical 
studies would also be very useful to verify that the bounce observed in homogeneous models 
remains also present in lattice LQC, as well as studying the matter energy density p. In 
homogeneous models, p is bounded above by the critical energy density p c = 0.41p P1 , and 
it would be interesting to see if the same holds in lattice LQC, possibly modulo some small 
fluctations. 

The canonical quantization in lattice LQC is more complicated than that of homogeneous 
models due to, among other things, the presence of the diffeomorphism constraint in addition 
to the scalar constraint. A potentially important point is that the scalar and diffeomorphism 
constraints are treated quite differently in lattice LQC. The simplest way to see this is 
by looking at the effective theory, where the diffeomorphism constraint is the same as in 
the classical theory, but the scalar constraint is significantly different. In other words, the 
classical diffeomorphism symmetry is not modified even when the space-time curvature nears 
the Planck scale. This might indicate that time-like and space-like diffeomorphisms are on 
a different footing in the quantum theory: the space-like diffeomorphism symmetry seems 
to continue to hold in quantum gravity, while the time-like diffeomorphism symmetry is 
modified when the curvature becomes large. It is not clear if this property is simply the 
artefact of a simple model or whether it which holds more generally, but it is certainly worth 
further thought. 

In this paper, we have seen how scalar perturbations can be studied in lattice LQC. How- 
ever, it is not as easy to put vector and tensor perturbations on a lattice. This is because the 
metric of a flat FLRW space-time with vector or tensor perturbations is necessarily nondiag- 
onal. Nondiagonal metrics are difficult to handle in LQC, even when they are homogeneous. 
This is because the holonomy of the Ashtekar connection around a square loop in this case 



is not an almost periodic function of the connection (see [14| for an explicit example of 
this in the Bianchi II space-time). If the holonomy around the square loop is not almost 
periodic in the connection, then it is not known how to build an operator corresponding to 
the field strength in LQC and then the Hamiltonian cannot be defined by the introduction 
of a nonlocal field strength operator as was done here. This problem can be avoided by 
introducing a nonlocal connection operator as in 14|, [l5j in which case a nondiagonal metric 
can be handled, thus making it possible to treat scalar, vector and tensor modes together 
in a gauge-invariant LQC treatment; we leave this for future work. 

A limitation of lattice LQC is that inverse triad effects cannot be included in the model 
— either in the quantum or the effective theories — as one of the necessary assumptions to 
work on the lattice in the first place was for each cell to be larger than the Planck volume and 
therefore also for the wavelengths of the modes of interest to remain greater than £p\, in which 



case inverse triad corrections are expected to be negligible 28]. Inverse triad corrections 
could play an important role, particularly in inflationary models where the wavelengths of 
the modes visible today were much smaller than £p\ at the bounce point in LQC but then 
lattice LQC is not an appropriate model to study such a scenario. Nonetheless, we expect 
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lattice LQC to be an excellent model for cosmological scenarios where the wavelengths of 
interest remain larger than the Planck length at all times, including the ekpyrotic model 
and various matter bounce scenarios. 

The main open question now is whether the quantum gravity effects studied here leave 
some signal in the classical regime in cosmology and, if they do, how they might be detected 
in for example the cosmic microwave background. This should be studied for the ekpyrotic, 
matter bounce and inflationary models (in which case inverse triad effects could also be 
important), as well as all other cosmological scenarios in order to precisely understand the 
observational consequences of quantum gravity effects in each of these cosmologies, and 
thereby test LQC. 
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